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THE CLASSICAL "SPHERE-OF-INFLUENCE" 

/ S' 

SUMMARY 


The classical definition of the "sphere-of- influence" is given and 
the equation implicitly defining its radius is derived. This equation 
was solved using a little known iteration technique attributed to 
Wegstein. Two analytic approximations are derived and the Fourier coef- 
ficients through 12th order for the exact result are displayed for the 
moon and each planet in the solar system. Several "spheres" about the 
moon computed according to these results are displayed graphically for 
the earth-moon system. 



I. INTRODUCTION 

It has been an increasingly popular and useful tool in the prelimi- 
nary design phases of interplanetary trajectories and even lunar probes 
to approximate the actual motion of a rocket vehicle by segments of 
trajectories found from the classical two-body equations. This means, 
for example, that the motion of a Mars probe might be approximated to a 
high degree of accuracy by three segments: an earth-centered conic 

segment, a sun-centered conic segment and a Mars-centered conic segment. 
The point at which the transition from one segment to the next occurs is 
somewhat arbitrary, but it has been customary to use an approximation to 
the "sphere-of-influence" as the criteria. This modern term, "sphere-of- 
influence," is a carry-over from the classical work of Laplace and 
Tisserand who used the term "activity-sphere." The surface to which 
both these terms refer is not truly spherical, but generically the terms 
are descriptive. It is not the intent here to discuss the effect of the 
"sphere-of-inf luence" on the accuracy of patched-conic techniques but 
merely to make available a reference on a much seen but little discussed 
term. 


II. DISCUSSION 


Considering n point masses, the equations of motion for the kth body 

are 


n , 



j=i 




X 


where K is the gravitational proportionality constant, 1 indicates that 
the dummy variable j cannot equal k, and the origin of the coordinate 
system is arbitrary except that it is not coincident with one of the 
masses. If the origin is translated to one of the masses, say the ith 
mass, the equations of motion can be written 


r k 


= -K(m. + mp 




I ft* I 3 





r>*| 3 


J 



( 1 ) 


where k^i, j^k, j/i and x indicates that the origin for the radius 
vectors is coincident with the ith mass. Notice for the classical two- 
body problem, equation (1) reduces to the well-known result 



KCm! + m 2 ) 



( 2 ) 


where m x is at the origin of the coordinate system. The right-hand side 
(R.H.S.) of (2) is the force per unit mass acting on m 2 and its magnitude 
will be called the two-body force. 

Now consider three bodies where the coordinate system in Figure 1 
is referenced. 


z' 



Figure 1 
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With the aid of (1), dropping the * notation, the equation of motion 
of m with respect to M becomes 


.. 

r* = -K(m + M) — - Km! 

|rf 


?-^| 3 l^l 3 


From Figure 1, it is seen that "p = r* - f* x and making the assumption that 
m « M, the last equation becomes 


r + Km 


1— >1 3 

-luni ^ 

y.|— »i 3 

1— » 1 

m 


1 PI 

l r xl 


The R.H.S. of (3) may be considered a perturbing force _^iue to which 
disturbs the motion of m about M. Calling this force F^., its magnitude is 


l?il • ?i - Km! 0r + ^ 

^ 1 ^ 1 


where p = |"p| and r x = |? 1 | . The ratio of the magnitude of the perturb- 
ing force to the magnitude of the two-body force for this case is 


m i (\ , 1 , 2p • ri 


Ri = r l + r* 


Using (1) again, the equation of motion for m with respect to m x is 


, K(m + mi) j?_ . km f 

i7?i a V- (-?!)i i-ni 3 - 


-K(m + m ± ) - KM 


+ ri 

-> , -» 1 3 

p + r il 
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Assuming m « m x and 


noticing that r* = "p + r^. 


this becomes 


■p + Km x 


JL 






( 5 ) 


Again, the R.H.S. of (5) may be considered a perturbing force due to M 
which disturbs the motion of m about m x . Calling this force F^ 2 , its 
magnitude is 



= km 1 : 


2r 


& 


1/2 


where r = |r*| and rx = l^il* The ratio of the magnitude of the perturb- 
ing force to the magnitude of the two-body force for this case is 


R 2 = p 2 


_M 

m l 


7 * 


1 

*1 ' 



( 6 ) 


The "sphere-of-inf luence" is now defined as the locus of points 
for which these two ratios are equal. Equating (4) and (6), evaluating 
the scalar products and rearranging yields 


p 4 = (mi/M) 2 r 4 


T 


p 4 + r 4 + 2p 
+ rj - 2r r 


2 

2" 

1 


2 

r 7 cos 9 

p cos 0 - 2r r 



( 7 ) 


where 0 is the angle between and This is the exact defining equa- 

tion for p, the radius of the "sphere-of-inf luence p is the only 
unknown in (7) since r = (r 2 4- p 2 + 2r x p cos 0) 1 / 2 , r x is supposed 
known, and 0 is a parameter. For m x < M, equation (7) represents a sur- 
face enclosing m 1 and p/r ± < 1. Equation (7) must be solved for p 
numerically using an iterative procedure. Wegstein*s iteration, dis- 
cussed in Appendix I is applicable and extremely rapid. However, 
numerical experience with a normalized form of (7) indicated that its 
denominator was not being evaluated accurately. The fundamental reason 
for this is due to the R.H.S. of (5) which requires differencing two 
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nearly equal quantities; always an inaccurate process unless special 
procedures are adopted. This computational problem is neatly circum- 
vented by applying an idea of Encke's which he employed in computing 
planetary orbits. The specific details of the application are dis- 
cussed in Appendix IX with the following result 


R 4 = (ma/MV 


(1 + R 2 + 2R cos 9) (1 + R 4 + 2R 2 cos eW 2 
[F 2 (l + R 2 + 2R cos 0) + 2F(R 2 + R cos 0) + R 2 ] 1 / 2 


= (mj/M ) 2 


(1 + Qi) (1 + R 4 + 2R 2 cos e) 1 / 2 
[F 2 (l + Qi) + 2F (Q i + R*) + R 2 ] 1 / 2 ’ 


(8a) 


where R = p/r x , Qi = R 2 + 2R cos 6 and F = F(Qi) is defined as 


F - F(Q X ) = - f Qi + X Qi + • • • + (-l) n ^ Q° + . 


n'. n'. 2 
n = 1,2, ..., if 0 I Qi < 1 


(8b) 


F « F(Q 2 ) = f iQa + ^ - ja Q| + & Q* + 


+ (-D 


(2n -5)1 1 n 

n'. (n - 3)'. 0 2n-4 ^2 


(8c) 


n = 3,4,5, ..., if -1 <Qi<0, 


where 


Q 2 - 


. Ql... 

1 + Qi * 


The Wegstein iterative technique was able to solve equations (8) to 
12 significant digits in at most 6 trials and in the majority instances 
3 or 4 trials for the moon and planets. This precision is not justified 
by the limited accuracy of the mass ratios involved, but assuming them 
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exact, the speed of convergence of the iteration scheme is as stated. 
Because the complexity of equations (8) requires an electronic computer 
or desk calculator for efficient evaluation, two simplified approxima- 
tions were developed in Appendices III and IV. These are given as (9) 
and (10) which explicitly define p. 


P = 


IL. 


(M/m!) 2 / 5 (1 + 3 cos 2 0) 1 / 1 °- y cos 0 * 


1 + 6 cos 2 0 
3 cos 2 0 


(9) 


P = 


l/lO _ 2/5 

2A fmj 

M 


,\ , s2/5 

1 + l |) (j) firi cos 0 " 50 cos 20 + 


= .91244 


2/5 


r l 


1 + .72995 


® 1 
M 


2/5 


cos 0 - .06 cos 20 + ... 


( 10 ) 

They are essentially variations of Tisserand's development. The maximum 
relative error incurred with (9) is less than 1.5 percent for the moon, 
negligible for the planets out to Mars, and less than .2 percent (usually 
much less) for the remainder. The maximum relative error using (10) is 
less than 51/2 percent for the moon and less than 31/2 percent for the 
planets. These two equations should be compared with the frequently 
seen expression for p 


p = (lUi/M ) 2 ^ 5 r !. 


( 11 ) 


Equation (10) suggests a Fourier expansion for p. Since p is an even 
function of 6, the expansion is of the form 


p(-fl) = a D + a x cos •a + a 2 cos 20 + . . . a n cos n0 + ..., 
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where 


a 

o 





o 


p(3) d$ 


and 



cos n9 d-6. 


for n = 1, 2, ... . 


The coefficients through 12th order were computed for the moon and 
planets and are displayed in Table I along with the largest relative 
error incurred in truncating the series at that coefficient. Trapezoidal 
integration was used for the calculation using the planetary data in 
Table II. A comparison of the actual coefficients with the first three 
coefficients yielded by (10) is made in Table III. 

For the earth-moon system, equations (8), (9), (10) and (11) were 
used to plot typical "spheres-of-influence" about the moon. These results 
are shown in Figure 2. Notice that these are cross-sectional shapes, and 
the three-dimensional surface is obtained by a revolution about r^_ or, in 
this case, the earth-moon- line. The maximum radius occurs at about 80.8° 
and is not on the earth-moon line as implied by Tisserand. Also shown 
for comparison is the gravity sphere defined as the locus of points for 
which the gravitational attraction of the moon equals the gravitational 
attraction of the earth. Its equation is 


P 


2 



( 12 ) 


which is the equation of a circle with radius 


^ 

ft -0 

and center displaced away from the moon's center by an amount 
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III. CONCLUSIONS 


Considering three bodies, two massy and one a particle of negligible 
mass, the "sphere-of-inf luence n has been defined as the locus of points 
for which the ratio found by taking the perturbing force of one massy 
body on the particle and dividing it by the two-body force between the 
particle and the other massy body is set equal to a similar ratio with 
the roles of the two massy bodies reversed. The resulting equation for 
the ,f sphere-of-inf luence" was conveniently and rapidly solved by using 
Wegstein’s iteration. Formulas for two analytic approximations and coef- 
ficients for Fourier expansions through 12 terms have been displayed. 
Which of these should be used depends on the accuracy required and the 
computational aids available. Finally, graphical display of some of 
these results indicates that, while the departure from sphericity of the 
"sphere-of-influence" is relatively small, It is nonetheless significant. 
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TABLE 
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TABLE II 


Solar System Parameters 


Planet 

ri(km)* 

( mi /M) 2/5 

(M/mi) 2 / 5 

M/m a* 

Mercury 

57.9 x 10 6 

.19313 x 10“ 2 

517.78 

6.1 x 10 s 

Venus 

108 x 10 s 

.57037 x 10 -2 

175.32 

4.07 x 10 5 

Earth 

150 x 10 6 

.61846 x 10" 2 

161.69 

3.3244 x 10 5 

Mars 

229 x 10 6 

.25352 x 10 -2 

394.44 

3.09 x 10 s 

Jupiter 

776 x 10 6 

.61937 x 10 _1 

16.145 

1.0474 x 10 3 

Saturn 

143 x 10 7 

.38226 x 10’ 1 

26.160 

3.5 x 10 3 

Uranus 

287 x 10 7 

.18064 x 10 _1 

55.358 

2.28 x 10 4 

Neptune 

450 x 10 7 

.19230 x 10" 1 

52.002 

1.95 x 10 4 

Pluto 

589 x 10 7 

.60585 x lO' 2 

165.05 

3.5 x 10 5 

Moon 

384.402 x 10 3 

.17212 

5.8099 

81.357 


These values taken from "Lunar Flight Handbook," NAS8-5031, 
Martin Co. , 1963. 

These values refer to earth-moon system. 
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TABLE III 


Comparison Between Analytic Fourier Coefficients from 
Equation (10) and Their Numerically Determined Values 


ORDER 



i 

2 

3 

SUN-MERCURY 

Approximate 

Exact 

102036. 

103235. 

143.852 

114.303 

-6122.15 

-6918.37 

SUN- VENUS 

Approximate 

Exact 

562072. 

568684. 

2340.19 

1861.44 

-33724.3 

-33111.0 

SUN-EARTH 

Approximate 

Exact 

846471. 

856430. 

3821.40 

3040.03 

-50788.2 

-57394.3 

SUN-MARS 

Approximate 

Exact 

529734. 

535962. 

980.326 

779.091 

-31784.1 

-35917.8 

SUN- JUPITER 
Approximate 
Exact 

43855400. 

44420300. 

1982780. 

1602230. 

-2631320. 

-2976600. 

SUN- SATURN 
Approximate 
Exact 

i 

49878700. 

50486000. 

1391820. 

1117080. 

-2992720. 

-3383550. 

SUN-URANUS 

Approximate 

Exact 

47305800. 

47866200. 

623787 . 
497853. 

-2838350. 

-3207880. 

SUN-NEPTUNE 

Approximate 

Exact 

78959600. 

79896000. 

1108380. 

884894. 

-4737570. 

-5354460. 

SUN- PLUTO 

Approximate 

Exact 

32560700. 

32943800. 

144000. 

114552. 

-1953640. 

-2207750. 

EARTH-MOON 

Approximate 

Exact 

60371.7 

61649.2 

7585.28 

6355.15 

-3622.30 

-4093.37 





FIG. 2. SPHERE OF INFLUENCE ABOUT MOON 









APPENDIX I 


A classical method for finding the root of a nonlinear equation is 
known as the method of iteration. If a root of F(x) = 0 is desired and 
F(x) can be solved for x in some form x = 0(x) , then the algorithm is 


*n = 0<*n-l>- 


(I-D 


A sufficient condition for convergence to a root is 


l0'(xn-l) I < 1 


(1-2) 


for each of the in some neighborhood of the root. The sequence defined 
by (1-1) is demonstrated graphically in Figure 3. The desired root, a, 
is obviously the intersection of the straight line and the curve 0(x). 

In this particular example, the sequence (1-1) is oscillatory, nonconvergent , 
principally because \0' (a) | > 1. Wegstein devised a modification of (1-1) 
that can force convergence even when (1-2) is violated. The modification 
consists of altering the successive values of x n . Denoting the altered 
values by xi, then 


*n = Xn + K (xn - i n _l>’ 


(1-3) 


where now Xjj = ^(Xjj,^). An estimate for K can be derived by referring to 
Figure 4. Ideally, K(x n - x n _^) = -BC, and since x n - x n _^ = AC, 


Now 


so that 


K = 


BC 

AC 


1 + K 


AC - BC = AB 
AC AC 


K = BC AC = BC 
1+K ' AC AB AB ' 
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From Figure 4, BC = BD, and by the mean value theorem, there exists a 
x n _^ < ^ < a such that 


0(a) - 0(x n _ 1 ) 
0'(i) = : 

a - x n _! 


BD 
AB * 


Therefore, 


K = BC 
1 + K ” AB 




or 


= 0'(S) 

1 - 0'(S) * 


(I 


0' (|) can be estimated as 


0'(i) 


0(x n -i> - 0(x n - 2 ) 

x n-l " x n-2 


x n - x n-l 

Vl • *n-2 


(I 


Substituting (1-5) into (1-4) into (1-3) yields 




Notice that 0 % (2) = -3/2. Actually, Figures 3 and 4 are plots of (1-7) 
and show that (1-7) unmodified diverges as expected. The classical 
iteration and the Wegstein procedure are compared in the following table 




1 

1.00000000000 

1.00000000000 

2 

2.87500000000 

2.87500000000 

3 

.029541015625 

.029541015625 

4 

1.74476338246 

2.99999677755 

5 

1.93922535110 

-.374989124219 

6 

2.00484760174 

3.00659122337 

7 

1.99991096179 

-.397294289403 

8 

1.99999987059 

3.00783875302 

9 

2.00000000000 

-.401524978457 

10 


3.00809184799 

11 


-.402383715948 

12 


3.00814387697 

13 


-.402560265342 

14 


3.00815460127 

15 


-.402596656742 

16 


3.00815681300 

17 


-.402604161995 

18 


3.00815726919 

19 


-.402605709998 

20 


3.00815736328 

21 


-.402606029282 

22 


3.00815738269 

23 


-.402606095173 

24 


3.00815738669 

25 


-.402606108735 

26 


3.00815738752 

27 


-.402606111558 

28 


3.00815738768 

29 


-.402606112111 

30 


3.00815738773 

31 


-,402606112257 

32 


3.00815738773 

33 


-.402606112257 
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APPENDIX II 


The numerical problem mentioned in the text comes about as a result 
of attempting to evaluate the magnitude of 



(II-l) 


where r = |r*| and r^. = (r^j . When r* and r\ are approximately equal in 
magnitude, the subtraction in (II-l) causes a loss of significant digits. 
This means that, if a certain number of significant digits are desired in 
the result, ^ and r x must be known with a greater number of significant 
digits (assuming r* and r^ are sufficiently close in magnitude). The 
requirement for extra precision can be by-passed using a modification of 
Encke's well known astronomical method. Substituting = r^ + p into 
(II-l) gives 


ri + P 
r 3 




or 


— * 




ri . 

73 p 


(II-2) 


The R.H.S. of (II-2) can be evaluated accurately if (r^/r) 3 - 1 can be 
found accurately. There are essentially two cases to consider: (1) 

r § rj. and (2) r < r l4 


Case 1 : Define Qi by 


1 + Qi = (r/r x ) 2 . 


(H-3) 
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(r x /r) 3 = (1 + Qo.)" 3/2 


= 1 - 2 ^ + “ 


q! + 


+ (-i) 


n (2n + 1) 


n'. n'. 


by the binomial expansion, valid for |q x | <1. Consequently, 


(r x /r) 3 - 1 = - 


4 Qi + 


Q^ + 


n » 1, 2, 


+ (-l) n ^ + 1 ^ > ' Q 
v / 2 n x 


n'. n'. 2 


= F(Q x ). 


Substituting (II-4) into (II-2) yields 


^3 |f(Q x ) + (l + F(Q X )) , 


where from (II-3) 


Ql = (r/ri) : 


. 1 - r? <?,+?>• -M3) - rf 


2i? • rN + f? * "p p 2 + 2p ri cos 0 

p = p 

*1 r i 


= R 2 + 2R cos 0 



Case 2 : Define Q 2 by 


1 + Q 2 = (rj./ r) 2 . 


(H-7) 


(r^r ) 3 = (1 + Q 2 ) 3 / 2 


= 1 +^Q £ 


+ 2 . Ifil. 2 . I . 1-2 1_ 


2 3-2 


by the binomial expansion, valid for |q 2 | < 1. Consequently, 


( r i/ r ) 3 - 1 = f ( Q 2 + 2 ^ " 3 Q2 2 3 + 4-^2 4 " 4 ^2* + ' * • 


c_^\ n — C^ n ~ 2iLi 1 — o 11 + 

( i} n'.(n - 3)'. 2 2 n -4 Q 2 + * * 


. , n = 3, 4, 5, 


= F(Q 2 ). 


Substitution of (II- 8 ) into (II-2) yields a form similar to (II-5) 


|f(Q 2 ) + (i + F(Q 2 )J 


where from (II-7) 


Qs = (ri/r) 2 . 1 = = sLz S^+Jl • (r 1 + a 


r (fl + p) • (*1 + p) 


(2^ • r^i + 1 ^ • ^ _ _ (P + 2P ri cos 6 ) 

-> . 0 -» -> r? + 4- 2 p r-, cos 0 

r i * r i + 2 P * r i + P ' P 1 1 
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(II - 10) 


Q2 ~ “ 


R + 2R cos 9 
1 + R 2 + 2R cos 


1 + Qi* 


Referencing (II-5) and (II-9) in the form 


jiFr x + (1 + F) pk 


the magnitude of the perturbing force of M disturbing the motion of m 
about mi is 


|f 2 | = . % = P |f ^ + (l + f) = P |f(?i + ■$) + c?| 


= ™ F 

* 1 L 


2 (r’i + p) • (?! + p*) + 2FC?! + ~P) ' "p + 


P ' P 


= ^ F 2 (r 2 + 2p r x cos 9 + p 2 ) + 2F(p 2 +pr x cos 6) + 

r i L 


F 2 (l + 2R cos 6 + R 2 ) + 2F(R 2 + R cos 0) + R^ 


™ F 2 (l + Qi) + F(Q X + R 2 ) + R 2 


so that R 2 becomes 


R 2 = — R 2 F 2 (l + Qi) + F(Qi + R 2 ) + R 2 
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Setting Ri = R 2 gives 


mi p f l . 1 . 2o • r*i 

-t r 2 + -5 + 


7* + 7* + '^3 73* ' = — R 2 

1 


M \p’ r " r~ / mi 


1/2 


ni /2 


F 2 (l + Qi) + F(Qi + R 2 ) + R 2 


, , p . p . „ _v /ri + o 4 + 2p 2 rf cos 0 N \ 

(m 1 /M)(r 1 + p + 2p r x cos 0) ^ u — p 4 ^4 — m J 


i /2 


— R 2 


F 2 (l + Qi) + F(Qi + R 2 ) + R 2 


l/; 


or 


R = (m x /M) 


2 (1 ± R 2 + 2R cos 9) (1 + R 4 + 2R 2 cos Q) 1 ^ 2 
[F 2 (l + Qi) + FCQi + R 2 ) + R 2 ] 1 / 2 


= (mx/M) 


2 (1 + Qi)(l + R 4 + 2R 2 cos eW 2 


[F 2 (l + Q x ) + 2F(Q X + R 2 ) + R 2 ] 


, 1/2 ’ 


which is equation (8a) of the text. To apply the Wegstein iteration, 
this is rewritten as 


(1 ± ± El ± isi COS 

[F 2 (l + Qi) + 2F(Q a + R 2 ) + 


eW 8 

R 2 ] 1 / 8 


The F(Qi) and F(Q 2 ) series are both alternating series; thus, the trunca- 
tion errors are less than the first term of the neglected remainder of 
the series. Simple calculations show that 22 terms of the F(Qi) series 
and 12 terms of the F(Q 2 ) series are sufficient to yield an error less 
than 1CT 8 if Q x and Q 2 are less than or equal to .4. Actually, for the 
planets, the maximum Q x and Q 2 are an order of magnitude smaller than 
this; therefore retaining the quoted number of terms yields a precision 
much greater than required. In using these series, is always calcu- 
lated, and F(Qi) is computed if Q x § 0 and F(Q 2 ) is computed if Q z < 0. 
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APPENDIX III 


Equations (9) and (10) are derived by making use of the fact that 
p/r-L <1. First, rewrite (7) as 



Substituting 


r 

r i 



and 


R = 


— yields 

r i 



R 4 jl ± (1 + R 2 + 2R cos 9) 2 - 2(1 + R cos 9) (1 + R 2 + 2R cos 9> l/2 
(1 + R 2 + 2R cos 0) 2 \/l + R 4 + 2R 2 cos 0 


Expanding the radical and applying the binomial theorem several times gives 


= R 4 (1 - 4R cos © + •..) v/r 2 (1 + 3 cos 2 0) + 4R 3 cos 0 + . .. 


= R 5 (l - 4R cos 0 


■f . . . ) *1(1 + 3 cos 2 0) ^1 + 


4R cos 0 
1 +3 cos 2 0 


+ . . . 


= R 5 (1 - 4R cos 0 + 




+ 3 cos 0 vl + 




4R cos 0 
1 + 3 cos 2 0 


+ . • . 
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- i 5 — / . 2R cos 0 

= R 5 (1 - 4R cos 0 + . . . ) v 1 +3 cos 2 9(1+ — + 

1+3 cos 2 0 


= R 5 \/l + 3 cos 2 6 ( 1 + 


2R cos 9 
1+3 cos 2 0 


- 4R cos 0 + . . . 


c- / 5 — - i 2R cos 0 - 4R cos 0 - 12R cos 3 0 

= R 5 vl + 3 cos 2 9 ( 1 + + • 

1+3 cos 2 0 


- I 5 — f ( 1+6 cos 2 0 ) \ 

= R 5 vl + 3 cos 2 0(1 - 2R cos 8 1 . . . ). 

\ 1+3 cos 2 0 / 


Thus , 


R 5 ■ Kr 





+ 3 cos 0 


( 1+6 cos 2 0 ) 

2R cos 0 f- . . 

1 + 3 cos 2 0 


-l 


or 


R = 


© 


2/5 


(1+3 cos 2 0 ) 


l/io 


( 1+6 cos 2 0 ) 

1 - 2R cos 0 h 

1 + 3 cos 2 9 


Therefore , 

, N 2/5 

■ ■© — — 

(1+3 cos 2 0 ) 

which, solved for R, gives 


1/10 


2 ( 1+6 cos 2 0 ) 

1 + 7 - R cos 8 1- 

5 (1+3 cos 2 0 ) 


R = 
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so that 


P = 


M_ 


2/5 


(1+3 cos 2 9) 


1/10 - (!) 


(cos 0) 


(1+6 cos 2 0) 
1+3 cos 2 0 


which yields equation (9) and differs from Tisserand's result only in 
that he expands the denominator in the next to last equation. 
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APPENDIX IV 


Equation (10) is derived as follows. In equation (8), the denominator 
through second order terms in R becomes R 2 + 3R 2 cos 2 9. Neglecting terms 
higher than first order elsewhere yields 


-6 


— ) (1 + 4R cos 6 + ...) - 

'/r 2 (1 + 3 cos 2 0) + .. 


or 

R 5 

so that 
R 


~ (1 + 4R cos 0 + ...) (1 +3 cos 2 0+ ...) 


, 2/5 


J (1 + 4R cos 0 + . . . ) 1/ ^ 5 (1 + 3 cos 2 0 + . . . ) 1 ^ 1 ° 


m. 


2/5 


R ~ ( ■ J (1 + j R cos 0 + . . . ) (1 + 3 cos 2 0 + . . . ) 


-1 /io 


Substituting cos 2 0 = ^ (1 + cos 29) gives 


An \ 2 / 5 

R = (jf) (1 + j R cos 0 + . . . ) ( (5/2) (1 + | cos 20) + . . . 


-l/io 


solving for R yields 




2/5 r 


(1+1 cos 20) 1/10 


, x 1 / 10 
2 \ 

pJ vs j v; 


(C\ /V 

\5 J V M 


2/5 


cos 0 + . • . 


-1 


29 





Substituting 



0.91244, j = 0.8, and ^ = .06 yields equation (10). 
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